HUPD0910, UT-10-07 

Effects of Large Threshold Corrections in Supersymmetric Type-I 

Seesaw Model 

Sin Kyu Kang^, Takuya Morozumi^, and Norimi Yokozaki^ 

^ School of Liberal Arts, Seoul National University of Technology, Seoul 139-473, Korea. 

^ Graduate School of Science, Hiroshima University, 

Higashi- Hiroshima, 739-8526, Japan. 
O, 
CN , ^ Department of Physics, University of Tokyo, Tokyo 113-0033, Japan. 

> . 

O . Abstract 

in 



in 

o 
o 



We investigate lepton flavor violating (LFV) radiative processes and the relic abundance of 
neutralino dark matter in supersymmetric type-I seesaw model. We carefully derive threshold 



y^ • corrections to the flavor off-diagonal elements of slepton mass matrix and up-type Higgs mass 

o,' squared and find that they can be large in the case of large -Bjy. We examine how the branching 

,^ , ratios of LFV radiative decays and the relic abundance of neutralino dark matter can be significantly 



affected by the large threshold corrections. Soft scalar mass squared parameter of up-type Higgs 
scalar is also affected by the threshold corrections. Since the higgsino mass depends on the mass 



(N 
> 

in . 

^T) . parameter for up-type Higgs, the LFV processes and the relic abundance of the neutralino dark 



matter are correlated with each other. We show that there are parameter regions where the 
predictions of the relic abundance of neutralino dark matter are consistent with WMAP observation 
and the branching ratios of LFV radiative decays are predicted to be testable in future experiments. 
We find that the masses of scalar supersymmetric particles are not necessarily small so that the 



X 

C^ ' branching ratios of LFV decays can be testable in future experiment, which is distinctive feature 

of this scenario. 

PACS numbers: 



I. INTRODUCTION 

Seesaw mechanism has been invented to explain the smallness of observed neutrino masses 
relative to those of quarks and charged leptons. Among several varieties of seesaw models, 
the simplest version is type-I seesaw model which requires the standard model gauge singlet 
right-handed(RH) Majorana neutrinos and the existence of a huge mass scale which can 
be near or similar to the scale of grand unification. Supersymmetric(SUSY) version of the 
type-I seesaw model not only inherits this feature but also stabilizes the electroweak scale 
without fine-tuning, and provides a natural candidate for a dark matter. The neutrino 
Dirac type Yukawa couplings in SUSY type-I seesaw model are flavor off-diagonal, which 
give rise to neutrino mixing observed by neutrino oscillation experiments. Thanks to those 
Yukawa couplings flavor off-diagonal elements of the slepton mass matrix are induced by 
radiative corrections even though we take slepton mass matrix to be diagonal at the high 
energy scale. These flavor off-diagonal slepton masses can enhance branching ratios for 
lepton flavor violating(LFV) decays such as r — ;■ fij, r — ;■ 67 and /x — )■ 67 compared to those 



in non SUSY seesaw models 



m 



In our previous work, we have observed that threshold corrections to Higgs bilinear terms 
mediated by RH sneutrino can affect the minimization condition for the Higgs potential and 
thus the fine-tuning may be reduced when the mass splitting of RH sneutrinos becomes 



large 



B. 



We have shown that such a large mass splitting can be originated from large value 



of the B-term for RH sneutrino, ^ 

C = --BNMRN*^ + h.c.. 

In addition, we have also shown that those threshold corrections can significantly affect the 
relic abundance of the neutralino dark matter in minimal supergravity scenario(mSUGRA). 
Thus some of the parameter space excluded by WMAP data in SUSY seesaw model without 
the threshold corrections can be consistent with WMAP data when we include the threshold 



^ The large value of -Bat, (i.e. Sjv ^ 'msoft ^ B{^ clectorweak scale)) can be obtained from the term of 
superpotential given by [9|, 



I 



with Mb, = X{X), Bn = {Fx) / {X). In the case that {X) - lO^^ GeV and (Fx) -- lO^i GeV^ which is 
the same order of SUSY breaking F-term, B^ ^ 100 TeV is obtained. 



corrections. However, we have not considered the flavor effects generated from the threshold 
corrections. 

In this paper, we show that the threshold corrections can give rise to sizable contributions 
to lepton flavor violating phenomena, and extensively discuss how the branching ratios of 
LFV radiative decays are correlated with the relic abundance of neutralino dark matter. 
We believe that this observation is new although there are several literatures where the 
contributions of -B-parameter and threshold corrections to lepton flavor violating phenomena 



have been studied 
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This paper is organized as follows. In section II, we derive the threshold corrections to 
Higgs and slepton mass squared parameters by using RGBs for corresponding parameters. 
We also calculate the flnite terms which are not included in the approach using RGE method. 
In section III, we investigate how the branching ratios of LFV decays and the relic abundance 
of neutralino dark matter are affected by the threshold corrections and correlated with each 
other. In section IV, we devote to the numerical calculation and present our results. The 
concluding remarks will follow in section V. One-loop RGEs for slepton and Higgs masses 
including threshold corrections are presented in appendix. 

II. THE THRESHOLD CORRECTIONS 

In SUSY type-I seesaw model, flavor off-diagonal elements of SU{2) slepton mass matrix 
can arise from radiative corrections mediated by RH neutrinos and sneutrinos even though 
the slepton mass matrix is taken to be flavor-diagonal at the high energy scale such as the 
GUT scale. These corrections are evaluated with the help of RGEs given in Ref. [2|] under the 
assumption that soft scalar masses and scalar trilinear couplings are universal at the high 
energy scale. When the mass splittings of RH sneutrinos are large, threshold corrections 
arisen from integrating out heavy sectors should be taken into account. As we studied in 
8j, large mass splitting of RH sneutrinos is originated from large value of i?-terms for RH 
sneutrinos. It turns out from our numerical estimation that those threshold corrections 
can dominate over the other radiative corrections to the flavor off-diagonal elements of 
SU{2) slepton mass matrix. It is also worthwhile to notice that both radiative corrections 
and threshold corrections to the flavor off-diagonal elements of slepton mass matrix can 
signiflcantly contribute to Higgs mass squared parametersjsl, |9|. 



Now, let us derive the threshold corrections to slepton and Higgs masses by integrating 
out one-loop RGEs in the case that the mass differences among three generations of RH 
sneutrinos are large. The superpotential in SUSY type-I seesaw model is given as 

W = Y,,^H^ ■ UE, + Y,,jU ■ H^N^ + ]^Mr,^N,N, - ^iH, ■ H^, (1) 

and soft SUSY breaking terms are written as 

Csoft = -ml..ElE, - ml^llL, - m%^^^NJN, - ml^H, - m^^HlH^ 
+ {mjj^jj^Hi ■ H2 + h.c.) - (^A^,ijLi ■ H2N* + h.c^ 

' -72 aV* aV* 



- [A,,,E, ■ U~e\^ + h.c.) - - {BIn:n; + h.c.) . (2) 

Redefining the chiral superfields, iVj, we can take Mr^ij to be diagonal as 

Mn^ij = MR^idij, (3) 

where M^^i are real, positive and assumed to be hierarchical, i.e. Mr^i <^ Mr^2 ^ Mr^^. 
Then, in the case of \Bfj\ <^ \M'^^- — M'^j], the mass eigenvalues of RH sneutrinos are 
approximately given as 

Ml, = m%^^^ + Ml,-\Bll (4) 

Here, the mass eigenstates, Nh^i and Ni^i are given by 

N: = ±c'^''^/\N,,-^N,,), (5) 

where 

$. = arg(i?^,). (6) 

The hierarchical structure of RH neutrino and sneutrino masses is presented in FigUl The 
flavor off-diagonal part of the mass terms for RH sneutrinos is given as 

^ = -In {<,AA^>j + ^IAAj + miuAA,j + ^kiA.Nh,) , (r) 



where 
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FIG. 1: Mass hierarchies of RH neutrino and sneutrino masses. Mji^i denote RH neutrino masses 
whereas Mh^i and M;^j are heavier and hghter sneutrino masses, respectively 

Let us calculate radiative corrections and threshold corrections to the slepton and Higgs 
mass squared parameters. We note that there exist two types of contributions, one is the 
logarithmic contributions which can be obtained from the RGEs and the other is finite 
contributions which are calculated in this paper by the diagrammatic approach. 

A. Leading log contributions 

The leading log contributions to slepton and Higgs mass squared parameters can be 
obtained from RGEs. Let us first consider RGEs for slepton masses and up-type Higgs mass 
without taking into account the mass threshold contributions of RH sneutrinos. The explicit 



expressions of the RGEs arisen from only the contributions due to the RH neutrinos and 



sneutrinos are given by for slepton masses 
'dQ 



QE 






2m^^(F;rJ).„ (9) 



and for up-type Higgs mass 

ciim 

{16n')Q^ = 2TT{Y:m%Yj) + 2Tr(F>|F;) + 2m^ Jr(r;Fj) 

+ 2TriA:Al). (10) 

We assume that soft SUSY breaking parameters are universal at the GUT scale, so that the 
following relations are valid, 



mi ,, = mt-,,, = ml6ij, m^^ = ml, A^^ij = AoY^^ij, Bf^ = B%MR^i5ij. (11) 



2 2 

If .. = m ,% . . 

L,i] N,ij 



By integrating the above RGEs (MTUl) between M/j^ < Q < Mqut, "we obtain flavor off- 
diagonal slepton masses and the radiative corrections to up-type Higgs masses, and their 
approximated expressions are given by 

6m}.. ^-y ^S^ Ul + 3ml) In ^, (12) 

k ' 

and 

^^^. ^ - L ^^^ (^0 + 3mo) In ^^. (13) 

Here, we have used the fact that when Q"^ < M'^j^ the chiral superfield N^ decouples, and 
thus Yy^ik is set to zero. 

Next, let us derive the corrections to the above relations (I12|13l) originated from RH 
sneutrino mass thresholds shown in FigJH We split a complex scalar Ni to two real scalars 



Nh^i and Ni^i as in Eq.Q and then re-derive RGEs by using the method given in Ref. 12 |. 



The complete results are presented in Appendix A. The terms of the RGEs which contribute 



to the derivation of the threshold corrections are written as 



djTfi- 

k 



dQ 



R,k 



7 _, [y,y,ik^iy,kj^ '' + ^u,ikyu,kj^ '"J Mpi^j, 



[9(0' - Ml,) - 0{Q' - Ml, 



+ E 



fc-i 



Y* \u-2 \yT j_\^ (y* r2 p-i^kY'^ -lV* u'i* pi^kY'^ \ 

'^ u,ik\'^kk\^u,kj ^ 2-^ Vv,ik'^k'k^ ^ v,kj ^ ^v^ik^kk'^ v,k'j) 



k'=l 



[e{Q' - Ml,) - e{Q' - M, 



Lk)\ ' 



(14) 



and 






dQ 



+ 2 ^ Re(y,:,,Aj,,e-*'=)M«,, {d{Q^ - Ml,) - e{Q' - Ml 



i,k 



E 



i,k 
X 



fc-1 






[^(Q2 - Ml,) - 0(g2 - M, 



(15) 



Note that the terms of RGEs given in Eq. (TT4l) and (TTSjl have been derived in the diagonal 
basis of Majorana mass matrix M^ by following the method appeared in 13|. Integrating 
the RGEs for the slepton masses given in Appendix A for the region of energy scale, Mf, < 
Q^ < M^,, we obtain the threshold corrections as follows: 



St,mh. ^ m^^iMl,) - mljMl,) 



L,ij 



4i^i,, + 42^L-, 



(16) 



where 



xk 2 _ ^ [oY* y'^ m2 MY* A'^ -^fcfc I ,4* y'^ ^fcfc \ (17\ 



'■R,k 



42^1,, 



1 ^^^ 1 

TTZi Z-^ \-^mkBkkyv,inXu,kj + -^mk^ kkyu,ikyu,mj) j^2 ' \^°) 



m=l 



The terms which contain both Ay and S^ in Eq. (ll7p agree with Ref. ll| whereas disagree 
with Ref. [lOJ . The authors of Ref. ll|] derived the threshold corrections by using the method 



of analytic continuation into superspace |l4J-|l 6| and the author of Ref.(lO| derived the cor- 
rections by diagrammatic calculation. In the derivation of the above equations, we have 
used the following relations, 



In 



Ml, 



I R2 I 

Mlk 



m 



+ 



N.kk 



Ml, 



1 \Bl,? 

2 Ml, 



+ 0(M-J), In 



Mju 
Mh 



I R2 I 

'M' 



+ 0{M-l). (19) 



Imposing universal condition for the soft SUSY breaking parameters at the GUT scale, 
S^f^m^j. . . is approximately written as 

1 



4^1,, 



87r2" 



Y*,,Yj;,^{ml + AoB%). 



(20) 



We note that since BfAi ^ j) are radiatively generated at one- loop level, Sfj^2''^} • ^^® two- 



L,ij 



loop contributions. To estimate how large it is, we solve RGE for Bf- given in Eq.( ]A5|) . and 
then obtain its value around at Q = Mfi,,3, 



Bl^iMns) = BIMr, 



mk 



(FjF;)^fe(l + <5„fc), Mgut 
m 



87r2 



M, 



R,3 



AoMii,k —^ m ■ 



4vr2 



M 



R,3 



(21) 



where m < k. We notice that the term proportional to B^ in the RH side of Eq. (!2T]) can 
be diagonalized by changing the basis of Mji into the diagonal basis at Q = M/j3, but the 
term proportional to Aq can not be simultaneously diagonalized as proven in appendix C. 



Then, the contributions to 5^/j2'^l ■ ■ ^^^ approximately given by 



B% Aq 

47r2 167r2 



/ ^ \^u,im\^u ^1/ Jm3^jy,3j + J^u,i3\^u ^u )3mJ^u,mj) 



,m=l 



In 



M 



GUT 



M 



R,3 



(22) 



which turns out to be much smaller than 5j\]^m| . ., so that we can neglect these contributions. 
Similarly, integrating Eq.( IT5|) for the range of the scale, M;\ < Q^ < M^f^, we obtain the 
threshold corrections to the Higgs mass fnj^^ given by. 



^Wh, 






kYuMV^O 



AoBl) 



(23) 



B. Scheme dependent finite terms 

There are also finite terms which are renormalization scheme dependent. We obtain them 
by calculating corresponding Feynman diagrams and estimate how large they are. We used 
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MS (DR) scheme to subtract divergences. It turns out that these terms are small compared 
to the logarithmic contributions. To derive the scheme dependent terms, we first consider 
the divergent diagrams which arise from RH neutrino-sneutrino Yukawa interactions. The 
Feynman diagrams relevant to our purpose are shown in Figs j2|3l If we insert more flavor 
off-diagonal slepton masses into those diagrams, they can lead to only finite corrections 



proportional to the mass parameters such as B^ 



m- 



.. and m^ . . with i ^ j^ which are at 



«J' "''L,ij """ "''N,ij 

most at two-loop level because the mass parameters are radiatively generated at one-loop 
level under the assumption that soft masses are universal at the GUT scale. 

The one-loop corrections to m?j .. are obtained by calculating the diagrams (a)-(c)(g)(h) 
in Figs. 2,3. The corrections are composed of two contributions arisen from scalar loops 
and fermion loops. 



'^L,ij 



L,ij 






(24) 



The scalar loop contributions are given by 






-j# 



k \ A* Y^ p 



i<J>t 



)Mr,, [/: 



2l"^22>K,fc) 



UmUM^ 



+ 2 ^ v,ik^u,kj^ 

+ ly:,kYj;kj[^himl,) + (l + \^^\yMl,)hiMl,) + (l + \fi\yMl,)hiMl,)] 



221 ^'"^l,kJ 
2 \j ( TL/r2 



+ ^:^l^,e^^~^'^ 



-^^liM'UKki M?,') + ^mlLk'kh{Ml,„ M?,) 



+ i 



(25) 



mikk'UKki Kk') + mUkMMtk. Ml',,) 

Z Y* Y'^ 

327r2 '*'^ ''^-'' 

where 77122 — l/^P + '^H2- '^^^ ^^^^ momentum dependent term of Eq. fl25|) is obtained by 

and by keeping the term which 



Ml7k 



expanding —iT,^- . . with respect to external momenta — 2- 
remains in the limit of large M^^k- The loop functions /i and I2 are given as, 

d'^k 1 



him') = Q 



4-d 



(27r)^ F - m' 



m 



2 r 



12(771^, 7n2) = Q 



4-d 



-1 + 1 - In 



d'^k 



771 

~Q' 



(26) 



(2tx)'^ k'^ — mf k'^ — m^ 



167r2 



6-^ + 1 + 



771^ 



771% 



771^ 



1 ^l 



7717, 



77lt, 



rrin 



771^ 



'"e^ 



(27) 



The fermion loop contributions are calculated with the fermion loop diagrams by keeping 
external momenta assumed to be small compared to M^^, so as to derive the contributions 



from wave function renormalization. The results are expressed as 






/i| 



-2 Y, y:,kYj,k, 1 + 17^ h{Mi,) - ^s 



fc=i 



Mlu 



L,if 



(28) 



where S? . . is momentum dependent part given as 



- ^., = t -^:^l,f^ (f-' + ^ - '■' ^) + ''('''Hi 



(29) 



The exphcit form of — iSj^ •• is shown in Appendix B. Then, the quadratic parts of effective 
Lagrangian for sleptons Lj can be written in MS (DR) scheme as follows, 






^^tr\Q'')diXwi3 



L,ij 



— m 



2(fcmax)/^2\ 



2(fcmax)/^2\ , r(A;max),^2 



("^max/^^^ 2 



m + Smf^T-'iQ') + 51,--^^ Z., + Shi'^L,, LIL,, (30) 



where the contributions of the loop diagrams mediated by heavy neutrino superfields, A^^i 
^fcmax ^^^ included and the parameters are given by 

"^niax -I 71 /f 2 



z 



V"^niaxj / A^2 



i^max. -I A/f '^ 



L,t3 ^ ^ 167r2 



fc=l 



(31) 



2_^\^v,ik^u,kj + '^u,ilXu,kj^^H2 + Yu,ikYu,kj''^N^kk) 



Ml, 



k=l 



Q' 



^' * T 9 max(M2^,M2^,) 

^ Z^ ^u,ik-^u,k'j"''N^kk'^'-^ n2 



k^k',1 



Q' 



^tr'^'u, 



IGtt 



-I /^'max "^niax 

V2 / ^ V i^.tfc i^.fc? ^ ^ u,ik^ u,kj"''H2 ^ / J '' v.ik'- v,k' y"' fq ^]j_ki 1 1 



k=\ 



k'=\ 



-1 "^max 

Xl«">axj 2. _ -*- >P 



C(*:max) 2_ 

^SI ^ L,ij - Y^ 



n^ 



k=l 



r)Y* v^ m^ -\- (V* A^ U^* -X- A* Y^ R^ 

'^^ v,ik^ u,kj"''N^kk ' \^ v,ik^u,kj'^kk ^ ^u,ik^ v^kj-^k 



M, 



R,k 



(32) 
(33) 
(34) 



Note that 830^^% ■ • and 6sifn'j . . are the scheme dependent and independent finite terms, 
respectively. One can normalize the kinetic term in canonical form by replacing Lj as, 

V^i^-^L, ^ L„. (35) 

With the replacement, the effective Lagrangian is written as. 



/^'vmax 

-'-eff 



d^L*d^Li 



m 



-^("^max J 



'("^maxj L/ 



\"^max j 



+ 5mf'-^-^ + d's^'-^ml + ^ifr^^ml),,!*!,. (36) 
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Now let us apply the matching condition at the scale Q = Mr^u^^^ to two effective theories, 
one of which contains the superfields, A^i ~ A^fcmax) ^'^^ the other contains only A^i ~ ^femax-i- 
The effective Lagrangian with the fcmax active superfields should be the same as the one with 
the fcmax — 1 active superfields at the matching scale; 

£f-)(g = Mk,,_) = C%Y-'\Q = Mn,k^J. (37) 



Then we can derive threshold corrections for the soft breaking mass for sleptons and for the 
Higgs particle. Applying the matching condition at Q^ = M"^/^^^^, we obtain the relation, 



m 



'^yf^max. J- J 



{Mn,k 



\"^max J-/ 



1 



m 



^{f^maix) 



X'^TaEix) L/ 



(M, 



R,kn 



+ C"(^max) 2 r("^niax J-) _ 2 



■'SD 



Here, we have used the relation. 



SI '"-L- 



V^^max j 



("'inax -Lj 



(38) 



Smf--\Mn,,^J = Smf---'\Mn,,^J. 



(39) 



Inserting the wave function renormalizations given as 



'("-max j 



Lnj 



A .-\^J_Y* Y^ rin— ^-1 



k=l 



327r 



1 



_ V^'max j 



Ml 






k=l 



Q' 



_ v'^i^a^/ 



\"'max J-j 



6 a + ■ 
^ 3271 



Ml. 



Y* . V 



Q' 



into Eq.(l38l). we finally obtain 



(40) 
(41) 

(42) 



m 



■^V'^max J-/ 

L,ij 



(M, 



R,kn 



."£r'(A'K.'=~) - 3^ 



"^ ^LiT" (^R'k,^^^)yu,nk^,J^u,k^^^j 



Y* Y^ ■m^^'^™'"'VA/fD 7 

1 



167r2 



I^,ikmiix ;^,fcmaxi -/V,fcniaxfcmax 









i^max -L 

■^"^ --T _2 



Z^ -^^-fc'i"^JV,fcn,axfc' 



fc' = l 



T^max J- 

-' y,ik"''N,kkma.J I'.'i-maxj 
fc = l 



+ 



3^2 [^K^k.. 



167r2 



y m- 

max ;^,Kmaxi -/V,fcniaxfemax 



+ fy*- y4'^ B^* + y4* ■ y^^ B^ 



T'max'i'niax 



Z^jl/Cniax ^'i^'niaxj /^niaxfi^n 



^i?,fc„ 



(43) 
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where we use Bfj = B^Mn^iSij and both scheme dependent and independent threshold 
corrections are included. If Bq is larger than the other soft breaking parameters, i.e., Bq > 
A^yjiriff, the scheme dependent terms in second bracket in Eg. (143 p are smaller than the 



scheme independent terms in third (final) bracket in Eg. (143 p . Therefore we can neglect the 
scheme dependent terms for the numerical analysis. We notice that the scheme independent 
terms are the same as those obtained from RG analysis given by Eg. (ll7p . 

Similarly, one can obtain the threshold corrections for bilinear part of the Higgs field H2. 
Calculating the diagrams (d)(e)(f)(i) in Figs. 2,3, we get one-loop corrections to mj^^ which 
are divided by two contributions as follows. The scalar loop contributions are given by 



^^H2 ~ 2 \^y,ik^i^,ki + '^u,ik^v,ki^R,k) 



Hmlu^Mi,) 



+ hijn 



LX.) 



+ i^(Y:^,AI^^^'^ + A;,,i;^,,e^*'=)M«,, h{mi ^^, M^,) - h{mi ^^, M/ 



+ 77n:..>^iJ2/i(m| J + h{Ml,) + h{Ml,)] 



, ml 



+ Yl KikYu,kjml,J2 {ml,. , , . .^.^ 

1 






|(*fc-*ft') 



+ ^<i,k'kHMik',Mi%)-^<i,kk'^^(Kk^^i,k') 



<,kkMMiu, M,%,) + mf,,,,/2(M,^„ Ml^,,] 

\ i P Y* Y^ 



327r 



2 u,ik Ujki' 



(44) 



The fermion loop contributions are given by 

3 



-zEi^ = -2J]y^,l^,,Jl(M^,)-^S^^, 



k=l 



where 



(45) 






E 

fc=i 



167r2 



Y Y n 

uAk u,kiy 



T 2 / --1 



1 Ml^ , 



(46) 



The explicit expression of the total contributions to m^^^ is shown in Eg. flB2p . Then the 
effective Lagrangian for Higgs {H2) field is given as. 



''-cff 



ZH,d^H;d>'H2 



- {rnf-\Q') + Smf-\Q') + Si'^^mj,^ + ^r^rn^,^) H2^ H2. (47) 
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where 



ZH^iQ'^ 



8mf--\Q'') 



fe=l 

^ r^^niax 



T n„ ^ifc 



Q" 



-11 



(48) 



167r2 

"-'max 



2_^\^u,ik^u,ki + ^u,ik^u,kj^^L,ji "*" ^u,ik^u,ki'^N ,kk> 



Ml, 



k=l 



Q' 



^ T 2 1 max(M|„M|;,,) 



fcT^A:',! 



Q^ 



(49) 



-1 "-max ",max 






fc=l 

n^niax 



k'=l 



1 I'm ax 

W2 Z^ 



IGvr^ 



"-max 

E 

fc=i 



fc=i 



oy^* "v^T 2 _|_ (Y* A'^ r2* -I- 4* V'^ r2 ^ 

^-' i^.ifc-' u,ki"'']\f kk ^ \^ v,ik-^v,ki'^kk ^ ^iy,ik^ u,ki'^kk/ 



mi . . log -^^ -mi .. log ^^ 



M, 



i?,A; 



/--■ 1 fi-7r2 '^.fcj Lji f^fc 



L.ij ° m' 



L,ii 



i¥=j 



IQn^ 



1 1 

mi . . — m,i .. 

L,jj L,ii 



(51) 



Using Eg. (1471). one can obtain the soft mass of up-type Higgs including scheme dependent 
and independent threshold corrections as follows, 
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(52) 



Here, similar to the case of slepton masses, the scheme dependent terms are smaller than 
the scheme independent terms. 
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FIG. 2: Contributions from scalar degrees of freedom 
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FIG. 3: Contributions from fermionic degrees of freedom 

III. LEPTON FLAVOR VIOLATING DECAYS AND THE RELIC ABUNDANCE 
OF NEUTRALINO DARK MATTER 

As we have seen in the previous section, the threshold corrections to slepton and Higgs 
masses can be so large that they dominate over the RGE running effects which are ap- 
proximately expressed in Eq.( IT2l) and ( TT3l) . In this paper, for our numerical calculation, 
we consider so called minimal supergravity scenario (mSUGRA) where gaugino masses, soft 
scalar masses and scalar trilinear couplings are universal at the GUT scale. In particular, we 
investigate the relic abundance of neutralino dark matter in the focus point region 
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19|, 



which is one of the regions where the relic density is consistent with WMAP observation. 
In the region, the lightest supersymmetric particle (LSP) is a mixture state of bino and hig- 
gsino and the annihilation cross section for LSP is enhanced due to the appropriate portion 
of the higgsino component. In the seesaw model without threshold corrections, the focus 
point region is significantly affected by the neutrino Yukawa sector and shift to the energy 
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scale 



'ar from electroweak scale when the right-handed Majorana masses are sufficiently 



large 20l422| . In SUSY type-I seesaw model based on mSUGRA, the radiative corrections 
drive the up-type Higgs mass squared to be more negative than that in MSSM if we do not 
include the threshold corrections. This feature is due to the presence of RH neutrino and 
sneutrino sectors. In this case, the higgsino mass |/i| is larger than that in MSSM, too. This 
in turn leads to larger relic abundance of neutralino dark matter compared to that in MSSM 
for fixed values of the soft scalar masses due to the small portion of the higgsino components 
in the LSP. However, if we include large threshold corrections mentioned above, the up-type 
Higgs mass squared is driven to be less negative which leads to smaller \fi\. Therefore the 
portion of higgsino in the lightest neutralino state becomes large, and the right amount 
of relic abundance of neutralino dark matter consistent with WMAP observation can be 
obtained even in the parameter space of MSSM excluded by WMAP data [8|. 

Since the threshold corrections mentioned above produce flavor off-diagonal slepton 
masses, they are new additional source of lepton flavor violating phenomena. Since large 
threshold corrections lead to large flavor off-diagonal slepton masses, we anticipate that the 
amplitude of lepton flavor violating processes such as r — )■ /i7, r — )■ 67 and /i — )■ 67 can 
be enhanced due to the new source of lepton flavor violation. Therefore there exists a ten- 
sion between the branching ratio of LFV decays and the relic abundance of neutralino dark 
matter. 

Higgsino mass parameter /i is determined by the minimization condition of the Higgs 
potential given as. 



1 2 I ,2 , rnjj^iml)- minimi) tan^/3 
-my = — u H TT^ 

2 ^ "^' tan2/3-l 



-L 2 I 12 , ""HiV'ZJ iiu^ ynuz) ^a.LL ,j 



In the limit of large tan /3, this condition can be written as 

1 



^m| ^ -|/ip - m]j^{m\). (54) 



The radiative corrections to Ta^^ are given as, 



Y* yT 
dm, ~ ^ 



k 



il + AoB% - {Al + 3ml) In ^' 



(55) 



where the first two terms of RH side are the threshold corrections and the last term is the 
RG running effect. We see from Eq.( !55l) that in the case of the large threshold corrections 
to m|^2 5 the first two contributions to Sm,"^^ dominate over the last one and thus Higgsino 
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mass parameter |/ip becomes small. This lowers the abundance of neutralino dark matter 
compared to that of MSSM. From the constraint from WMAP observation, we can obtain 
the allowed parameter space. 

Now, let us consider the radiative LFV decays and investigate how they are related with 
the relic abundance of neutralino dark matter. Branching ratios of the radiative LFV decays 



a 



are approximately given by [3| 

BR{k ^ Ijl) ~ 7^^^ tan^ /3, (56) 

where rris is average SUSY scalar mass. For numerical calculations, we have used the com- 
plete formulas which are given by Eq.(51) in Ref.p]. As we can see from Eq. (13^ . the 
branching ratios of the LFV decays are proportional to \m'j ..p. Including the threshold 
corrections, m| . . can be approximately written as 



Y* yT 

rrif .. ~ > '■ TT-^ 

L,i3 A^ 87^2 

k 



ml + AoB% - {Al + 3m^) In 



Mgut 
11, 



R,k 



(57) 



Here, the first two terms of RH side are the threshold corrections given in Eq. (l20p . The 
last term comes from the RG running effects which is obtained by integrating out RGEs for 
m| . .. Because the relative sign is opposite, both contributions cancel each other when their 
contributions are comparable. As the B^ becomes large, the first two terms gets dominant 

?° 



over the last term. It is worthwhile to notice that the size of B'^ is limited by the constraints 



from the experiments. 



IV. NUMERICAL CALCULATION 



In this section, we present our numerical calculations of the relic abundance of the neu 
tralino dark matter and the branching ratio of the LFV decays. For our numerical calcula- 



tions, we use the micrOMEGAs package |23l.l24| for the relic abundance and the SuSpect |25l| 



for RGE running of soft SUSY breaking parameters with appropriate modification. Then 
we compare our numerical results with the experimental ones. We use the following up- 



per bound for the branching ratio of the LFV decays : BR(r — )■ fi'j) < 4.4 x 10 ® |26| . 



BR(r -> 67) < 3.3 x IQ-^ [26| and BR(/i -^ 67) < 1.2 x 10"" |27|. 



The WMAP observation leads to the relic abundance of the cold dark matter, VLcdm^ 
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2 



given by 28 1 



^CDuh — 0.111_q;oi5 



{2a) 



For numerical calculations, we need to parameterize the neutrino Dirac Yukawa couplings 
in terms of light neutrino masses and mixing. From the superpotential, the Lagrangian for 
the neutrino sector is given as, 

^ = -\yJ,kMh ■ H2 - Ir^.kl't ■ ^2iVr - p^Mn^kK' + h.c. (58) 

After integrating out the heavy Majorana neutrinos, one obtains the dimension five operator. 

^eS = +o(^j ■ H2)Y^,jkYJ ^u,kSi ■ -^^2), 

The mass terms for the left-handed light neutrinos are given by 

C = --uIM^^ul + h.c. , 



(59) 



(60) 



where the left-handed light neutrinos ly are presented in the flavor eigenstate, z/^ 
(z/e, z/^, z/,-), and Mi, denotes the neutrino mass matrix written as 



i-vT /rT0\2 



The neutrino mass matrix M,, can be diagonalized by the unitary matrix U,y as follows: 



(61) 



U'^MyUy = diag{mi,i, 1711,2, rn„3) = M^d- 



(62) 



The unitary matrix Uu can be identified to Umns, defined by 

/ 

Umns = 



C12C13 



S12C13 



Sue 



-iS 



\ 



\ 



~Si2C23 — Ci2S23Si3e* C12C23 — ^123235136* S23C13 

■S12S23 ~ Ci2C23Si3e* — C12S23 — Si2C23Si3e* C23C13 



X 



diag (. 



,i«i/2 ia2/2 



!)■ 



(63) 



where 6,ai,a2 are CP-violating phases, Cij = cos(%) and Sij = sin(6'jj) with mixing angles 
6ij. The current experimental values and bounds for the mixing angles are 



sin2(2^i2) = 0.87 ±0.03, sin2(2^23) > 0.92, sm^{29r3) < 0.19. 



(64) 
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Since we do not consider the CP-violation, the phases are set to zero. In our numerical 
calculations, we take 612 = 0.6, sin023 = cos 6*23 = l/V^ and sin ^13 = 0. It is well known 
that the neutrino Dirac Yukawa matrix can be written in terms of Umns and the diagonal 
forms of the mass matrices as follows |3|] ; 

yif = (^)^v^i?0i4^^1,^5, (65) 

where R is a complex orthogonal matrix. We assume that R = 1 for the sake of simplicity 
in our numerical calculations. 

The combination (Y*Yj)ij can be written as 

(y;rj),, = .^U,^^km,kMn,kUlkr (66) 

Assuming that neutrino masses are subject to normal hierarchy, we can take 



{m^i, m^2,m^3) = (0, y Am^^, y Am^^ + Aml^), (67) 

where Am^i and Am32 are mass difference of neutrinos, and their experimental results are 
Aml^ = (7.57 ±0.20) x lO^^eV^ and Aml^ = (2.43 ±0.13) x lO^^eVl We also assume that 
Mr^i < Mji^2 < M^,3- In this case, 

iY:yJh - j^U.,r3m.sMn,3Ul,3r (68) 

(-"2/ 

Thanks to sin6'i3 = 0, {YJYj)is <^ {Y*Yj)23. As we can see from eq.f l57|) . this implies that 
"^i 13 ^ "^123' Therefore, we anticipate that BR(r — > 67) ^ BR(r — ;■ ^^7) despite the 
current experimental bounds are of same order. So, we do not present BR(r — )■ ej) in this 
work. 

FigJUshows how both the relic relic abundance and the branching ratios of the LFV decays 
simultaneously depend on the parameter B^ when light neutrino masses are hierarchical. We 
take universal soft scalar mass rrio to be ITeV. Here, note that lightest SUSY particle is the 
lightest neutralino. We also take tan/3 = 5(10) in the upper (lower) panels. The other input 
values of the parameters we take are presented in the caption of FigJH The left (right) two 
panels in Fig^BR^r — > fij) {BR{jj, -^ 67)) vs. the relic abundance as a function of B% is 
shown. In each panels, green solid and red dotted curves present the branching ratio of LFV 
and the relic abundance of the lightest neutralino, respectively. The gray dotted (upper) and 
blue solid (lower) horizontal lines show the current upper bound on the branching ratio from 
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experiments and the relic abundance of the dark matter obtained from WMAP, respectively. 
When B^ is so large that the prediction of ^x^"^ ^^ ^° ^^^ observed abundance of the dark 
matter, BR{t — )■ /i7) and BR^jj, — )■ 67) are predicted to be quite large. In this case, the 
prediction of BR{t — )■ yU7) almost reaches to the experimental bound. Although there will 
be a chance to probe the LFV decays in the case of large value of S^ in future experiments, 
the size of B% is limited by constraint coming from the relic abundance of neutralino dark 
matter. In FigJSl we also show how the branching ratios of LFV depend on the universal 
soft scalar mass ttiq. Here, there exist the valleys where the branching ratios are extremely 
suppressed, which are occurred due to cancellation between RG running effects and the 
threshold corrections. This is expected from Eq.f l57p . since the sign of the contribution from 
RG running effect to the soft SUSY breaking terms for slepton is opposite to the threshold 
correction. Therefore as mo increases, the large threshold correction from B% is cancelled 
and branching ratios for LFV become small. As can be seen from Fig. [5l the dotted (dashed) 
line starts from the point corresponding to itiq = 700(1400) GeV below which electroweak 
symmetry is not radiatively broken due to too large corrections to rrijj^ . In the case of large 
B% such as 400 ~ 600 TeV, the branching ratio for r — )■ /i7 and fi ^ e'y become so sizable 
that the LFV decay could be detected even when rriQ is larger than ITeV. 
In the case of inverted hierarchy, three neutrino masses are given as 



{m^3,m^i,m^2) = {m^s, yml^ + ^mj^ - Am^^, yml^ + Amj^), (69) 

and 777,1,3 < '^i/i < 'mu2- The threshold corrections to Smfj^ is proportional to Ty{Y*YJ). 



Ti{Y*Yj ) is approximately given as 



Tr(i;%^) ^ [(|f/.,i2r+|f/.,22r+|t/.,32r)m,2M^,2 

+ (|f/.,23|' + K,33\^)m,3MR,s'\ I {Hlf . (70) 

In the case of m,y3/m^2 ^ ^i?,2/jW^R^3, the first term in Eq. (!70|) is dominant. The ratio of 
Tt{Y*YJ) in the inverted hierarchy to that in the normal hierarchy is 

Tr(y;if )/Tr(y;yJ)^^^ ~ Mn,2/Mn,3. (71) 

Therefore threshold corrections to frijj^ are much smaller than in the normal hierarchy case. 
This leads to rather larger abundance of the neutralino dark matter. On the other hand, 
the term (Y^^YJ)^^ is written as 

(y;if )^^ ^ U^,22Ul,,m^,2Mn,2/ {Hlf . (72) 
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This term is obviously larger than that in the normal hierarchy. In particular, for 171^3 ~ 0, 

(>;%"),,/(n%^)rr~io. (73) 

Therefore we naively expect that BR{^ — )■ 67) becomes larger by 2 order of magnitude than 
that in the normal hierarchy for the same value of B%. As a result, it becomes difficult 
to satisfy the constraint from the current bounds on BR{fi — )■ 67) and the relic abundance 
simultaneously. 

When the neutrino masses are almost degenerate, i.e. m^i ~ m^2 ~ f^us, (Xu'^J)2i 
is larger than that in the normal hierarchy, which makes the prediction of BR{fi — )■ 67) 
in this case substantially enhanced compared to the normal hierarchical case, and thus 
the constraint of BR{fi — )■ 67) becomes more severe when we consider the constraint of 
the relic abundance simultaneously. We present the branching ratios of the LFV decay 
yU — )■ 67 and the relic abundance of the lightest neutralino as functions of B% in FigJHl The 
different figures correspond to different values of tan/3. As tan/3 increases, larger value of 
B^ is preferred to accommodate both BR{n — )■ 67) and the relic abundance of dark matter 
candidate. 

V. CONCLUSION 

We have investigated the LFV radiative decays and the relic abundance of the neutralino 
dark matter in the SUSY seesaw model. We have carefully derived the threshold corrections 
to the flavor off-diagonal elements of slepton mass matrix and up-type Higgs mass squared 
and found that they can be so large in the case of large B^ that the branching ratios of the 
LFV decays and the relic abundance of neutralino dark matter can be signiflcantly affected. 
Our numerical results show that there are parameter regions where the prediction of the 
relic abundance of neutralino dark matter is consistent with WMAP observation and the 



branching ratios of 
future experiments 



FV radiative decays can be enhanced so as for them to be probed in 



29 



30| . Although the origin of such large B% is unclear, if such large 
B-term exists, the branching ratios of LFV decays are signiflcantly enhanced even when mo 
is not small. Therefore, the masses of scalar supersymmetric particles are not necessarily 
small so that the branching ratios of LFV decays can be testable in future experiment, which 
is distinctive feature of this scenario. 
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FIG. 4: The branching ratios of LFV decays and the reUc abundance of the hghtest neutrahno are 
shown as functions of B^ with hierarchical neutrino mass case. Two left panels show BR{t — )• /U7) 
vs. the relic abundance and two right panels show BR{fi — t- 67) vs. the relic abundance. In each 
panel, green solid and red dotted curves represent the corresponding branching ratio and the relic 
abundance of the lightest neutralino, respectively. The gray dotted (upper) and blue solid (lower) 
straight lines represent the upper bound on the branching ratio given by experiments and the 
abundance of the dark matter obtained from WMAP, respectively. We take {Mji^i, Mjj^2; -^ij.s) = 
(10^'^, 10^^, 10^^)GeV and tan/3 = 5 for upper two panels and tan/3 = 10 for lower two panels. We 
also take mo = ITeV, Aq = 300GeV, mi/2 = SOOGeV and /i > 0. 
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FIG. 5: The branching ratios of LFV decays are shown as functions of universal soft scalar mass 
niQ for hierarchical neutrino mass case. The left panel shows BR[t — )• ^^y) and the right panel 
shows BR{p, — )• 67). The gray dotted lines correspond to the experimental upper bound on the 
branching fractions. Green solid, red dotted and blue broken lines are calculated with B^ = 
(100, 400, 600)TeV, respectively. We choose mi/2 = 300GeV,^o = 300GeV,tan/3 = 10. 
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Appendix A: Renormalization Group Equation 

In this appendix, we present one-loop renormahzation group equations for SU{2) slepton 
masses and Higgs mass squared parameters including the threshold effects. The one-loop 
RGE for the parameter Bf, is also presented. Here, we have omitted the same contributions 
as in MSSM. The RGEs for SU{2) slepton masses are given by 






E 



fc-i 



Y* 2 yT ,\^(Y* ™2 yT , y* 2 yT 

^u,ik"''N^kk-^i',kj ^ Z-^y-^i',ik'"''N,k'k-^i^,kj ^ ^ u,ik"''N^kk' -^ i^,k', 

k'=l 
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where m^g — "^^2 "'" l**^!^' "^^^ RGEs for Higgs mass squared parameters are given by 

fc-i 



[16.')q"< 



dQ 



E 

i,k 






X [e{Q' - Ml,) + ^(Q^ - Ml)] 
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i,k 

i,k 
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i,k 



+ 2 ^ ReiY*„Al,,e-^''^)Mn,k [9iQ' - Ml,) - 9iQ' - M; 



i,k 



+ E 



i,k 



fe-1 
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i,k 



(A3) 



where m^^^ = ||Up + m|^^. We note that the RH-side of Eq.( ]A3[l corresponds to the 
contribution from wave function renormahzation. 



;i6.^)Q^ 



i,k 

+ Y.f,Y*„A^„[e{Q' - Ml,) + 0{Q' - Ml)] 

i,k 



i.k 



The RGE for Bf, is given as 



[lQn')Q^ = 4 
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Ai,ky:,kjMR,j + mr,yi^,a,,u, 
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The RG equation for Mr can be found,for instance, in [4| 



fl 



;i6vr^)g^ = 2 {yjy:)^^ m«,,, + 2m^,, {yjy:) 



dQ 



jk ' 



(A6) 



Appendix B: Diagram calculation 



By summing the self-energies from both fermion and scalar loop diagrams, the one loop 
contribution to the slepton mass squared is given by. 
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The contribution to the up type Higgs {H2) mass squared is given by, 
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T^ '^^ v,ik^ v,ki"''N,kk ^ \^ iy,ik^u,ki kk ' ^i/,ik^ v,ki'^kk) 



+ /_^\^u,ik^u,k>iBkk'Bk>k' + '^u,ik>'^u,ki k'kBk'k'] 



Mr, 



k<k' 



Mlk^i 



+ T7-^n%n> r ' + 1 - In 



^L, 



IGtt 



Q' 



(B2) 
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Appendix C: Approximate solutions of the renormalization group equations for Mpt 
and B'^. 

As stated below Eq. fl2ip . although the renormalization group running may induce the 
large flavor off-diagonal contribution to Bf- (i 7^ j) at a lower mass scale, by switching 
the basis of Majorana mass matrix Mr to the diagonal basis, we can keep Bf, in the basis 
almost diagonal because off-diagonal elements BfAi 7^ j) are doubly suppressed by a factor 
of -g^ and one loop suppressed factor. The effect of the small off-diagonal elements in Bf- 
on slepton soft breaking term turns out to be smaller than the leading threshold corrections 
given in Eq.(l20l). 

In this appendix, we first show that the large radiatively generated off-diagonal elements 
of Bfj corresponding to the second term of Eq. (12T]) are indeed rotated away in the diagonal 
basis for M^. To show this, the renormalization group equations for B^ and M^ are solved in 
perturbative way, i.e., we use the approximation so that in the RH side of the renormalization 
group equations, all the couplings A^, , Y^, and mass Mr are scale independent constants 
defined at GUT scale where the initial conditions for renormalization group equations are 
imposed. We also show that off-diagonal elements of the third term of Eq. (pT]) remain even 
after the rotation and are numerically small compared with the leading diagonal elements. 

The solutions for Eq. (]A5|) and Eq. (]A6l) aX Q = M^3 with the boundary conditions in 
Eq.(!TTl). are given as, 

Mffi,(M«,3) = Mr^Aj - 2{H,,Mr,, + H,,MR,i)t'\ (CI) 

4(M^,3) = Bo[MR,i6i, - 2{HijMR^, + H,,MR,)t''] - 4Ao{H,,Mr,, + Hj.Mr,)^^ 

= iBo + 2Ao)MRiMR,3)i,-2AoMR,Aj, (C2) 

where t°^ = i^log^^^ and H = Y^Y;. Since the first term of B'^ in Eq.Q is pro- 
portional to the running mass matrix of the heavy Majorana neutrinos, it is also diagonal 
in the diagonal basis for Mr. As for the second term proportional to Aq, it is changed 
into the non-diagonal one. To derive the unitary matrix O diagonalizing Mr{Q = Mr-^) 



27 



approximately, we first write the mass matrix a.t Q = Mr,^ in the matrix form as, 



Mr{Q = Mh,3) 

-2H,2t'^^{Mn,i + MR,2) 



Mr2{Mr^^) -2i723t°'(Mfi,2 + Mh,3; 

-2if23t°'(Mfi,2 + Mh,3) Mra(MK3) 



\ 



(C3) 



where we have used the property Hij = Hji since CP is assumed to be invariant. The 
diagonal elements at Q = M^j 3 are given by 



03 \ 



MRi{Q = Mr,:,) = MR,i{l - AHui! 



One can find the matrix O with which the mass matrix Eq. (lC3p is diagonalized as, 



3J, 



0(M^,3)M^(Q = Mh,3)0' {Mr,,) = D{Q = Mr 



(C4) 



(C5) 



where D is the diagonal matrix. The rotation given above corresponds to changing the basis 

N'^ — > O^N'^. In the new basis, B^ is given by. 



Bl^,^ ^ {OB'O' 



{{Bo + 2Ao)Di{Q = Mr,,) - 2AoMR,,]5i^ 



2AoJ2i0^k{MR,k-MR,s)0^ 



jk) 



(C6) 



fc=i 



From Eq. (lC6p . the off diagonal elements of -B^ew ^^^ given as, 

2 

Bl^ij = -2^ Yl 0^k{MR,k - MR,3)0,k {t ^ 
fe=i 



J) 



(C7) 



The diagonal elements are dominated by the term proportional to BqMr^i, 

3 



Bl^^^ = BoD,{Mr,,) + 2Ao(A -J^OIMr,,). 



(C8) 



fc=i 



To write the off-diagonal elements of B^^^ explicitly, we introduce the parametrization for 
the orthogonal matrix O^ as. 



/ 



O' 



^N„N 
''13''12 



''13''12 



'13 



\ 



_^N „N„N _ „N „N _„N „N ^N ,„N„N „N „N 
''23''13'"12 '-23''12 ''23''13'='12 "r '-23'-12 '-13''23 

'-23"'13'-12 "T *'23''12 '-23''13''12 "'23'-12 '-13'-23 



(C9) 
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where s^, = sin^j^ and cf^ = cos9f^. One can write -Bnewii (^ ¥" J) using the angles as 



-2Ao{444(M^,i - M^,2) + 44(44 + 444)(Mij,2 - m«,3)}, 



R2 

-'-'new, 23 



r2 

-°new,12 



iV2, 



-2Ao{44((Mh,1 - Mr,2) - 4 (^R,l - ^R,3)) 



^.N I N N N , 



-,N\2\ 



+ 4(4424(1 + (3-)^) - C-3.f3 COS 2^- )(M^,2 " M«,3)}. (CIO) 



The angles 9j^ can be determined by the diagonalization Eq. (1C5l) . When Mr^s ^ Mr^2, ^/?,i, 
one can determine s^, s^ from the equation, 



/ 



-2Hrst''Mn,3 



\ 



-2H23t''MR, 







y-2H,st^'Mn,3 -2ff23t°3Mfi,3 MnsiMns) J 
Then we find that D^ ~ Mr^^, and s(^ and s^ are given as 

-2Hrst'', 



( s^ \ 

■'IS 

''13'='23 
. ^.N N 



D. 



■^13 

'-13'^23 

y^l3^2Zj 



.(Cll) 






'13 

4 ^ -2H,;t'\ (C12) 

where we ignore the corrections of the order of 0{Tf^) (i = 1,2). The determination of s^ 



^Mn- 



is more involved. It is determined by diagonalizing the following 2x2 matrix which can be 
obtained after the largest eigenvalue state is decoupled from the 3x3 matrix in Eq.f lC3p . 



^12 


■^12 


„7V 
■J 12 


4 


Di 


o\ 





D2 



Mr,i - Mn,3sj^ 
-2ifi2t°=^Mp 2 - Mp 3S13S23 



^N „Af 



-2i7i2t"'MR,2 - Mfl,3Si3S23 \ / CI2 SI2 



Mr,2 - Mfi,3si3 



'■^12 '-12, 



s]^ is approximately given by, 



42 ^ -2i7i2t°' - 



M, 



M, 



R,3 N N 
■'13'' 23' 



R,2 



(C13) 



(CM) 



where the following conditions are assumed to be satisfied, 

M^,2 > {s^3?Mr,3, (4)'^K,3- (C15) 

Using the formulae given in Eq. (1C12[) . one can write the dominant terms for the first two 
equations in Eq. flClOp . 



B, 



new, 13 



2AoMp3(4s 



N N 
12''13 



^N N^ 
■'12"' 23^) 



N „N I „N „N^ 



54^,23 = 2AoM^,3(44 + 4^ 



13;- 



(C16) 
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When s{2 ^ 1, they are simphfied as, 



-Snew,m3 " SAoM/j^sS^g 



-4AoM^,3i^m3t"', (m = 1, 2). (C17) 



Eq. (]C17p shows that in the diagonal basis of Mr-, the off diagonal elements of 5^3 is given 



by the third term of Eq. (l2Tl) and it is small compared with the large diagonal element 

Next we show that the variation of the threshold corrections Eg. (1201) due to the change 
of the basis is two loop effect and thus negligibly small. When we change the basis as, 

[NrY -^ O'^iV^^ (C18) 

Yy in Eq.( 12U]) should be replaced by 

n ^nO(M^,fcf. (C19) 

Then the threshold correction is replaced by the following equation, 

^t>L, = -^,iY:0^).u{OYj),,{ml + AoBl). (C20) 

Now let us examine how large 5^-^rn?^^.. for the case /c = 3 could be after changing the basis 
as follows, 

3 2 ^, ^0 + ^0-^JV ry* yT , y* (yT ^N , yT ^N -. , ry* ^N ,y* ^N -.yT N| 

(C21) 

Since s^3 = —{YjY*)m3,^t^^ , ijn = 1,2), the variation 5^j^m|.. is the two loop order and 
thus small correction. 

For the completeness, we check whether the conditions Eq. (lC15ll and s^ <^ 1 which 



lead to Eq. flU17p are satisfied. For numerical estimation. We take Mqut = lO^^(GeV), 
Mh,3 = lOi^(GeV), Mr^2 = 10i2(GeV),and Mr^i = 10^°(GeV). For Hij, one may use the 
parametrization Eq. (l65!) . 

H^, = iYjY:% = ^^^MR,MR,{Rm,R%. (C22) 

When R is a real orthogonal matrix, 

l-^*il < n .''2 o'^ ^ (l"^^3 -myi\ + \my2 -mi,i\). (C23) 

v^ sm p 
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Using tan/3 = 5 and v = 246(GeV), we obtain 



4 



\Hu\ < 1.0 X 10-^ 
\Hi3\ < 1.0 X 10-^ 
\H23\ < 1.0 X 10"^ (C24) 

Since 2t°^ ~ 0.06, we obtain 

|s^| < 6 X 10-^ Isfgl < 6 X 10-^ (C25) 

Then Eq. OClSp is satisfied. We also note that s{^ is as small as 

Isfal < 6 X lO"*' + 3.6 x 10"^ ~ 1.0 x 10~^ (C26) 
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